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This is a short non-technical introduction to applications of the Quantum Field Theory 
methods to graphene. We derive the Dirac model from the tight binding model and 
describe calculations of the polarization operator (conductivity). Later on, we use this 
quantity to describe the Quantum Hall Effect, light absorption by graphene, the Faraday 
effect, and the Casimir interaction. 
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1. Introduction 

Graphene , whic h is a one-atom thick layer of carbon atoms has many exceptional 
propertieJU^El making it one of the most interesting topics in condensed matter 
physics. The principle feature of graphene is that the quasi-particle excitations sat- 
isfy the Dirac equation, where the speed of light c is replaced by the so-called Fermi 
velocity vf — c/300. Therefore, the quantum field theory methods are very useful 
in the physics of graphene. By applying these methods, one can explain anomalous 
Hall Effect in graphene, the universal optical absorption rate, the Faraday effect, 
and predict the Casimir interaction of graphene, and do much more. 

The purpose of this article is to give a non-technical introduction to and a short 
overview of the use of quantum field theory in graphene. We start in the next Section 
with a derivation of the Dirac model from the tight binding model and a discussion 
of possible generalizations of the former. Quantum filed theory calculations in the 
Dirac model are presented in Sec. [3] at the example of polarization operator. This 
operator is then used in Sec. 14.11 to explain the anomalous Hall conductivity of 
graphene, which is proportional to n + 1/2 with an integer n. Next, in Sec. 14.21 we 
derive the universal absorption rate of air ~ 2.3% for optical frequencies. Sec 
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Fig. 1. The honeycomb lattice of graphene. 



is devoted to explanations of the measured giant Faraday effect in graphene and to 
further study of the polarization rotation. Sec. 14.41 contains a survey of calculations 
of the Casimir interaction of graphene. 

2. The Dirac model 

The Dirac model for quasi-paxticles in graphene was elaborated in full around 
198zJU5]_ twenty years before actual discovery of graphene. However, its basic prop- 
erties, like the linearity of the spectrum, eta, were well known and widely used much 
earlier due to the 1947 paper by Wallace^. The purpose of most of the works of 
the time was to describe graphite rather than graphene. For more details on the 
development and formulation of the Dirac model we refer the reader to a recent 
review ^ . 

In graphene, the carbon atoms form a honeycomb lattice (see Fig. [T]) with two 
triangular sublattices A and B. The lattice spacing is d — 1.42A. The nearest 
neighbors of an atom from the sublattice A belong to the sublattice B, and vice 
versa. By adding to the location of an atom in the sublattice A any of the three 
vectors 

*i=d(-l,0), da = d(l/2,V3/2), 5 3 = d(l/2, - V5/2) (1) 

one arrives to the location of one of the nearest neighbors in the sublattice B. 

In the tight binding model only the interaction between electrons belonging to 
the nearest neighbors is taken into account, so that the Hamiltonian reads 

3 

H = -tJ2 ^2(a^r a )b(r a + ^ + &t(rQ + ^XO) ' ( 2 ) 

aeA j=l 

where t is the so-called hopping parameter, and the operators (V , a, b\ b are creation 
and annihilation operators of electrons in the sublattices A and B, respectively. We 
adopt the units h = c = 1. These operators satisfy usual anticommutation relations. 
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Let us represent the wave function through a Fourier transform. 

|V) = (Vu(k) E ^-^(ra) + ^(k) £ e^&tfo)) |0) (3) 

The sublattice A is generated by shifts along the vectors S 2 — <5i and 5s — S\. There- 
fore, two momenta ki and k 2 are equivalent if (111(62 — 5\) + 712(63 — 8\)) ■ (ki — k 2 ) € 
27rZ for all integer n\ and n 2 . Representatives of the equivalence classes can be taken 
in a compact region in the momentum space - the Brillouin zone, which is a hexagon 
with the corners at the points V\ = 2w/3d (l, l/\/3), v 2 = 2n/3d (l, — 1/\/3) , 
t>3 = 2n/3d (0, — 2/V3) , «4 = —v\, t>5 = — v 2 and «6 = — V3. Opposite sides of this 
hexagon are identified, and, in particular, the corners v±, U3 and t> 5 are equivalent 
between themselves, as well as v 2 , V4 and v 6 are. 
One can calculate 

3 3 
Hty) = -tfa E E e!kr " flt ( r «) + ^ E E e^fo)) |0) (4) 

3=1 a£A j=l /3eB 

so that the stationary Schrodinger equation H\tp) — E\ip) becomes the matrix 
equation 

U-T)(£)-*(£)- < 5 > 

Clearly, the eigenvalues read 

£ = ±i|X|. (6) 

The spectrum is symmetric, with positive and negative parts meeting at the points 
where X = 0. Solutions to this condition are easy to find, and they coincide with the 
corners of the Brillouin zone. As we have discussed above, only two of the corners 
are independent. Let us take K± = ±v & = +v 3 as two independent solutions, called 
the Dirac points, describing two independent ground states. 

The next step is to expand the wave functions around these Dirac points, 
^ B (q) = iPa,b(K± + q). We suppose here that |q| is small compared to l/d <~ 
IKeV. Then, one obtains the Hamiltonians 

3td ( iqi ± q 2 \ , s. 
H ± = — . = v F (-a 2 qi ± crig 2 ) , (7) 

2 V -igi ±q 2 y 

where is the Fermi velocity, <jj are the standard Pauli matrices. By substituting 
in v F — (3td)/2 the value of d given above and t — 2.8eV one obtains a number 
which is slightly below the commonly accepted value vf — 1/(300). By introducing 
a four-component spinor -0 = (tp^,ipg,tp^,ip^), we can unify H± in a single Dirac 
Hamiltonian H — —ivF-j°j a d a , a = 1,2, where we replaced the momenta iq by 
partial derivatives, and 

^-0?£). M'M)- < 8 > 
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These 4x4 gamma matrices are taken in a reducible representation which is a direct 
sum of two inequivalent 2x2 representations. 

Each of these two-component representations for graphene quasiparticles' wave- 
function is somewhat similar to the spinor description of electrons in QED3 + i. 
However, in the case of graphene this pseudospin index refers to the sublattice 
degree of freedom rather than the real spin of the electrons. The whole effect of the 
real spin, which by now did not appear in presented Dirac model, is just in doubling 
the number of spinor components, so that we have 8-component spinors in graphene 
(N = 4 species of two-component fermions) . 

We conclude, that the tight binding model is equivalent below the energies l/d 
to a quasi- relativistic Dirac model, where the speed of light is replaced by vp. One 
should keep in mind that the tight binding model itself is an approximation. One 
can improve this approximation by considering, for example, couplings between 
atoms in the same lattice A or B, that is a next-nearest neighbor coupling. The 
corresponding energy is estimated to be about 0.1 eV, much smaller than the nearest 
neighbor coupling t considered above. Other possible couplings are summarized in 
Ref. |U It is expected that with suitable modifications the Dirac model is valid at 
least until the energies of ~ 2 eV. 

Under various circumstances, it may make sense to consider the following mod- 
ifications of the Dirac model. 

• One can add an interaction with the electromagnetic field. To preserve 
gauge invariance, this interaction is introduced by replacing the usual par- 
tial derivatives by gauge-covariant ones: d — > d + ieA. The electromagnetic 
potential is not confined to the graphene surface, but rather propagates in 
the ambient 3 + 1 dimensional space. This field may be an external magnetic 
field, a classical electromagnetic radiation, or quantized fluctuations. 

• Quasiparticles may have a mass. This mass is usually very small, but the 
introduction of a mass parameter may be convenient on theoretical grounds, 
e.g., to perform the Pauli-Villars regularization. 

• A more important mass-like parameter is the chemical potential /i, which 
describes the quasiparticle density. This parameter can be easily varied in 
experiments by applying a gate potential to graphene samples. Without a 
gate potential, fi is usually very small in a suspended graphene, but can be 
significant for epitaxial graphene due to interaction with the substrate. 

• Impurities in graphene are described by adding a phenomenological pa- 
rameter r which reminds an imaginary mass and enters the propagator of 
quasiparticles exactly as e in the Feynman prescription of contour integra- 
tion. 

• Finally, most experiments with graphene are done at rather high tempera- 
tures, which can be taken into account by usual rules of real or imaginary 
time thermal field theory. 
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3. Polarization operator 



Let us proceed with quantum field theory calculations based on the Dirac model. 
From the Dirac Hamiltonian one can derive the action 

S = J d 3 x$lfiip , Jj) = ij j {dj + ieAj) + ... (9) 

where the dots denote any of additional terms described at the end of the previous 
section. Tilde over 7 means that the space-components are rescaled 

7 =7°, 7 1 ' 2 = vpj 1 ' 2 . (10) 

Quantized fermions give rise to an effective action for external electromagnetic field 
(given by a sum of one-loop diagrams). To the second order in A, it reads 

S cS (A) = A w^^w A = 1 J ^P_ M . p)W i ip)Mp) ^ (n) 

where IP is the polarization operator. 

As an example, let us consider the polarization operator for a single massive 
two-component fcrmion at zero temperature, zero chemical potential, and without 
external magnetic field. Simple calculations give 



rrmn _ rn 
~ „2 'lj 



V? (12) 



with e J being the Levi-Civita totally antisymmetric tensor normalized according 
to e 012 = 1, T]J = diag(l, vf, vf), P m = VnP n - a = e 2 /(47r) ~ 1/137 is the fine struc- 
ture constant. The tensor structure of n m ™ is uniquely defined by quasi-relativistic 
and gauge invariances up to two functions, and 4>, which read 

2mp — (p 2 + 4m 2 )arctanh(p/2TO) 

V(P) = 7p > (I 3 ) 

2p 

= 2marctanh(p/2m) _ % 
P 

here p = +-*JpjP' ' , and we assume m > 0. From Eq. (|12j) one might get an impression 
that the polarization tensor is singular at the vp — > limit, or, that IP 70 is greatly 
enhanced by the smallness of vf- A more careful analysis shows that it is not true, 
the polarization tensor remains finite in the vp — > limit provided the frequency is 
non- vanishing, po 7^ 0. 

For more complicated external conditions, like a magnetic field, or for a non- 
zero chemical potential, the quasi-Lorentz invariance is broken, and the polarization 
tensor has a more complicated form than Eq. (|12j) . 

Particular importance of W k is due to the fact that this tensor defines how 
the electromagnetic field propagates through graphene. The full quadratic action 
for electromagnetic field is — \ fd 4 xF^ v + S c s, where the effective action (JTTJ) is 
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confined to the surface of graphene, which we place at a; 3 =0. The equations of 
motion following from this action contain a singular term 

0„F"" + 5{x 3 )W p A p = 0. (15) 

We extended LT to a 4 x 4 matrix with U 3 ^ = IP 3 = 0. The equations (fT5")) describe 
a free propagation of the electromagnetic field outside the surface x 3 — subject 
to the matching conditions 

Afj,\ x 3 =+Q = A /J | x 3 = _o, 

{d 3 A^\ x 3 =+0 - (0 3 ^)Ua = _ o = n^^U 3=0 (16) 

on that surface. 

The polarization tensor and, more generally, Feynman diagrams involving 2 + 1 
dimensional fermions were considered in a large number of papers. Still in 1980's 
the function iff (for Vp = 1) was calculated in Ref. while the pseudotens or part 
was discussed about the same time in the context of the parity anomal y^O^ i 
A renormalization group approach to theories with vp ^ 1 was suggested in 
Ref. Q21 In t his century , extensive calculations were done by the Kiev group and 
collaborator J 13 | 14 | 15 | 16 j The formulas ([j^ (jjgj and ^ are cons i stent with that 

calculations. 



4. Physical effects 

We proceed with considering some quantum field theory effects in graphene. Specif- 
ically, we will be interested in applications of the polarization tensor IP- 7 . We shall 
see, that this quantity indeed defines important and interesting physics. 

4.1. Quantum Hall Effect 

First of all, we note that the polarization tensor can be interpreted in terms of 
the conductivity of graphene. Indeed, variation of the effective action (fTTj) with 
respect to Ak produces an expectation value of the electric current in graphene, i.e. 
jk ro n fei A;. On the other hand, in the temporal gauge, A — 0, the electric field E a 
with the frequency w is related to the vector potential by E a = iu)A a . By definition, 
the conductivity is a matrix relating j and E. In this way, one arrives at the relation 

<r ab = — , 0,6 = 1,2. (17) 

To study the conductivity tensor at zero frequency (the dc conductivity) one 
puts the graphene sample in a constant magnetic field perpendicular to its surface 
and measures the anti-diagonal conductivity as a function of the chemical potential 
/i. This is prec isely the set-up for Hall experiments. For a one-layer graphene it was 
observec S^D^ that the off-diagonal (Hall) conductivity is quantized according to 
the law 



Cl2 



(n+i), n = 0,l,2,... (18) 
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i.e., the conductivity is proportional to half- integer numbers. This particular type 
of the Hall Effect is called anomalous, or unconventional integer, or half-integer 
Quantum Hall Effect. Starting with the Dirac model, the behavior ([18]) was pre- 
dicted in Ref. [19] by using the Feynman diagram approach described above, and in 
Ref. [20] from numerical simulations. This wonderful agreement between theory and 
experiment was the first confirmation of existence of the Dirac quasi-particles in 
graphene. 

It is interesting to note, that the half-integer Quantum Hall Effect is observed 
in the mono-layer graphene only. The double-layer graphene, for example, exhibits 
conventional Integer Quantum Hall Effect. This does not however imply that the 
Dirac model is not applicable to double-layer graphenes. The later case may be 
recovered if one carefully considers phases of determinant of the Dirac operator^ES. 

4.2. Absorption of light 

Another physical effect defined by the polarization operator is the absorption of 
light by a suspended monolayer graphene. In the simplest set-up one can neglect 
the chemical potential, suppose that there is no external magnetic field and put 
temperature to zero, T = 0. Consequently, one can use the polarization operator 
(TT2"]) . where, because of N = 4 generations of fermions in graphene, the scalar part 
is multiplied by N, while the pseudo-scalar part <fi cancells out. The cancellation 
occurs due to the form of gamma-matices ([5]), containing two inequivalent repre- 
sentations related by the parity transformations. In other words, we need to make 
a substitution 

f^f N = N^, 0^0. (19) 

Let us consider a plane wave with the frequency to propagating along the x 3 -axis 
from x 3 = — oo with the initial polarization parallel to x , which is being reflected 
by and transmitted through the graphene sample 

g— iut J 3 + (fxx e x + T X yey)e 3 , x < (20) 

I (j'xx^x txy^y)^ 3 > X ^> 

where b XjV are unit vectors in the direction a; 1 ' 2 . The mass-shell condition (free 
Maxwell equations away of the graphene sample) implies = u). For such waves 
the matching conditions (1161) simplify, 

A I — A I 

/1 °lx 3 =+0 ^ a \x 3 = -0 

(d z A a ) x3=+Q - (d z A a ) x3= _ = a^ N (k)S b a , (21) 



where we used (T12j) and ([T9|) . The transmission coefficients can be easily found, see 
e.g. Ref. (23] 
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The intensity of transmitted light thus reads 

I=M 2 = l + ^^ + 0(a 2 ). (23) 
to 

One can reformulate this result in terms of conductivity by noting that a^> = iuia xx , 
as follows from §T2§ and (fTD) . 

At large frequencies, uj 3> 2m, we have ^ ~ —muj /4, yielding aim ^ n /lo ~ — ont. 
The same conclusion is als o y alid if the polarization tensor is calculated with more 
general external conditions^! Therefore, we confirm the prediction of Refs. [551 HEl 
[57J US] made on somewhat different theoretical grounds of universal absorption rate 
of air ~ 2.3%, that was confirmed by the experimentPS. 

Clearly, this uniform absorption rate is much larger than one would expect from 
a one-atom thick layer. 



4.3. The Faraday effect 

The Faraday effect reminds very much the Hall effect at "non-zero frequencies" , and 
this analogy was used in Ref. [SUlto conjecture that the former should be common for 
Hall systemfl The set-up is, therefore, very similar: a graphene sample subject to a 
constant magnetic field perpendicular to its surface. Instead of the Hall conductivity, 
we shall be interested in the rotation of the polarization plane of a light beam 
passing through the surface of graphene. The frequency of photons shall be kept 
as a variable parameter, as well as the chemical potential. It can be shown that it 
is usually sufficient to consider the zero-temperature case only, but impurities are 
essential. 

By solving again the matching conditions dT6l) for plane wave ([20]) , but now with 
a polarization operator calculated in presence of constant magnetic field, one finds 
that the angle 9 of polarization rotation and the intensity I of transmitted light are 
given by 

9 = -^^ + 0(a 2 ), l=l-Rea xx + 0(a 2 ) : (24) 

where we used (|17p to express the n components through diagonal and Hall con- 
ductivities of gra phen e. 

An experimental] made recently demonstrated a "giant" Faraday rotation angle 
of about O.lrad peaked at low frequencies for a magnetic field of about 7 Tesla. A 
theoretical study^ shows a good agreement between this experiment and the Dirac 
model. Besides, the Dirac model predicts other effects, like step-function like behav- 
ior of the rotation angle and the peaks at higher freque nci es. Another interesting 
theoretical observation is that although the data of Ref. EH can be nicely fitted by 
the Drudc formula for conductivity, this Drude-like behavior cannot be uniformly 
extended for all frequencies. 



a In Ref. 23 the Faraday rotation was related to possible non-compensation of parity-odd parts of 
the polarization tensor between various generations of fcrmions. 



November 21, 2011 1:18 WSPC/INSTRUCTION FILE qfextll'3 



Quantum Field Theory in Graphene 9 

4.4. The Casimir effect 

The Casimir effectP^, which is one of the main topics of this Workshop, is sometimes 
understood as any manifestation of the zero point energy. We consider the Casimir 
effect in a stricter sense, as an interaction of two uncharged well-separated bodies 
due to quantum fluctuations of the electromagnetic vacuum. In the framework of the 
present review, let us take a suspended graphene sample separated by the distance 
a from a parallel plane ideal conductor. Under these conditions, one can neglect m, 
/i, and r, but the temperature will be non-zero, in general. 

The lowest-order diagram which gives the Casimir free energy ii^l 

J x ~ ; (25) 

where the photon propagator satisfies conductor boundary conditions on the sur- 
face x 3 — a. (We use the free energy since this is the relevant quantity at finite 
temperature). In this expression, one of the boundaries (conductor) is taken into 
account exactly, while the other one (graphene) - perturbatively, at the first order of 
a. A better approximation may be obtained by considering a closed loop of a prop- 
agator satisfying both conductor boundary conditions at x 3 = a an d th e matching 
conditions (|T6|) at x 3 = 0. This boils down to the use of the Lifshiti^ formula 

* = T E j f£ M(l ~ e" 2 ^r»(l ~ e-^rSrSr)] , (26) 



n— — oo 



„(1.2) 



where p\\ = lj% + p 2 , and u) n = 2nnT are the Matsubara frequencies, 
the reflection coefficients for the TE and TM modes at each of the two surfaces. For 

(2) (2) 

the second surface, which is an ideal conductor, we have rvj^j = 1, — — 1. The 
reflection coefficients for graphene are calculated similarly to the case considered 
in Sec. 14.21 though at non-zero temperature and arbitrary tangential momenta the 
calculations are more complicated. 

It is an interesting exercise to check that at the a 1 order the Lifshitz formula 
reproduces the two- loop diagram (|25|) . At zero temperature both approaches give 
consistent result of the order of 2.7% (for the Lifshitz formula) of the Casimir inter- 
action between two ideal metald^In the view of Sec. 14.21 this is not very surprising. 
Some unexpected features appear at non-zero temperature^. The perturbativc re- 
sult (|25|) rapidly becomes unreasonably large for growing T signalling that we have 
left the perturbative region. Roughly speaking, this effect is caused by a compe- 
tition between two small parameters, a and vf- The non-perturbative free energy 
([26)) also grows with increase of the dimensionless parameter aT, and, at the large 
T asymptotics we have 

TC(3) 
16na 2 

which is just a half of the interaction between two ideal metals in the same regime, 
or is the same value as for non-ideal metals described by the so-called Drude model! 



^oo^-^, (27) 
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Thus, the Casimir interaction of graphene at high temperature is extremely strong. 
This agrees qualitatively with Ref. [37|where the Casimir interaction of two graphene 
samples was considered. 

We conclude this Section with some references. Other papers which study the 
Casimir effect for graphene are Refs. [38l [39l [40] S ome ea rlier calculations used the 
hydrodynamic model for the electrons in gr aphene™. Since this model does not 
reproduce the linear dispersion law characteristic for graphene, this line of research 
was abandoned. More details on the presents status of Casimir effect in graphene 
can be found in Ref. S3J 

5. Conclusions 

The main message of this paper is that quantum field theory calculations based on 
the Dirac model of quasiparticles are extremely effective in describing the physics 
of graphene. One of the reasons for this effectiveness is the equivalence between 
the tight binding model and the Dirac model for small momenta. It is interesting to 
note that a single one-loop diagram of the polarization tensor considered in Sec.[3]is 
responsible for many physical phenomena, such as the Hall and Faraday effects and 
the uniform light absorption rate (where all experiments are in a good agreement 
with theory), and the Casimir interaction of graphene (where no experiment has 
been done so far). All the effects discussed above are very strong, much stronger 
than one would expect from a one-atom thick layers. 

Since parameters of the Dirac model may differ considerably from sample to 
sample of graphene, it makes sense to perform various types of experiments with 
the same samples. E.g., one can combine optical measurements with Casimir exper- 
iments. 

Some topics were not considered here, though they definitely deserve being men- 
tioned. One of such topics is the graphene nanoribbons. Before calculating the po- 
larization tensor, one should define boundary conditions which are compatible with 
quantum field theory. Such an analysis was performed in Ref. 1441 Another extremely 
intersting topic is the topological effects in graphene (see Refs. l45l |46| for a review), 
which includes the Jackiw-Pi models, applications of the index theorem, curvature 
effects, etc. This list of missing points is not exhaustive. There is much more in the 
area of applications of Quantum Field Theory to graphene. 
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